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Approximate analytical solutions describing the skyrmions given by rational map ansaetze are 
obtained. At large baryon numbers these solutions are similar to the domain wall, or to spherical 
bubbles with energy and baryon number density concentrated at their boundary. Rigorous upper 
bound is obtained for the masses of RM multiskyrmions which is remarkably close to known masses, 
especially at large B. The main properties of bubbles of matter are obtained for arbitrary number 
of flavors. 

1. Among many known soliton models used in different fields of physics the chiral 
soliton approach, starting with few basic concepts and ingredients incorporated in the model 
lagrangian provides, probably, the most realistic description of baryons and baryonic 
systems. The latter appear within this approach as quantized solitonic solutions of equa- 
tions of motion, characterized by the so called winding number or topological charge which 
is identified with baryon number B. Numerical studies have shown that the chiral field con- 
figurations of lowest energy possess different topological properties - the shape of the mass 
and B- number distribution - for different values of B. It is a sphere for B = 1 hedgehog , 
a torus for B = 2, tetrahedron for 5 = 3, cube for B = 4, and higher polyhedrons for greater 
baryon numbers. The symmetries of various configurations for B up to 22 and their masses 
have been determined in |3| (the references to earlier original papers can be found in |3], |] ) . 
These configurations have one-shell structure and for B > 6 all of them, except two cases, 
are formed from 12 pentagons and 2B - 14 hexagons, in carbon chemistry similar structures 
are known as fullerenes [Qc. 

The so called rational map (RM) ansatz, proposed for SU(2) skyrmions in || and 
widely used now allows to simplify the problem of finding the configurations of lowest energy. 
For the RM ansatz the minimization of the skyrmions energy functional proceeds in two 
steps: at first step the map from S 2 — ► S 2 is minimized for SU(2) model (for SU(N) model it 
is a map from S 2 — ► CF^ -1 , ||), and, second, the energy functional depending on skyrmion 
profile as a function of distance from center of skyrmion is minimized. As will be shown here, 
just the second step can be done analytically with quite good accuracy. Many important 
properties of RM multiskyrmions can be studied in this way, and some of them do not depend 
on result of the first step. This allows to make certain conclusions for arbitrary large B and 
for any number of flavors N F = N independently on presence of numerical calculations. 

So far, the chiral soliton models have been considered as a special class of models. 
Their connection with other soliton models would be of interest, and this is also an issue of 
present paper. 

2. Here we consider the SU(N) multiskyrmions given by rational map ansaetze; de- 
tailed comparison of analytical results with numerical calculations is made in SU(2) model 
and also in SU(3) variant using the projector ansatz ||. In SU(2) model the chiral fields 
are functions of profile / and unit vector ri, according to definition of the unitary ma- 
trix U e SU(2) U — Cf + isfrir. For the ansatz based on rational maps the profile / de- 
pends only on variable r, and components of vector n - on angular variables 9, <j>. ni = 
(2ReR)/(l + \R\ 2 ), n 2 = (2ImR)/{l + \R\ 2 ), n 3 = (1 - |i?| 2 )/(l + |i?| 2 ), where R is a rational function 
of variable z = tg(9/2)exp(i4>) defining the map of degree M from S 2 — * S 2 . 
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where f2 is a spherical angle. For B = 1 hedgehog Af = X = 1. AT = B for configurations of 
lowest energy. 

The classical mass of skyrmion for RM ansatz in universal units 37r 2 -F OT /e is H |6|: 

M = I- / (^vr 2 f 2 + 2i? S 2 (/' 2 + 1) +z4 V' (2) 
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r measured in units 2/(F w e), where we inserted the coefficient A N = 2(N — 1)/N for symmetry 
group SU(N) || to have a possibility to consider models with arbitrary number of flavors 
N = Nf - essentially nonembeddings of SU(2) in SU(N). The expressions for the quantities 
Af.l for projector ansatz in SU(N) model are presented in §. To find the minimal energy 
configuration at fixed Af — B one minimizes 1, and then finds the profile f(r) by minimizing 
energy (2). The inequality takes place: 1 > B 2 || ||. Direct numerical calculations have 
shown, and the analytical treatment here supports, that at large B and, hence, large 1 mul- 
tiskyrmion looks like a spherical bubble with profile equal to / = tt inside and / = outside. 
The energy and 5-number density of this configuration is concentrated at its boundary, sim- 
ilar to the domain walls system considered in |7j in connection with cosmological problems. 
Denote 4> = cosf, then the energy (2) can be presented as 
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with 4> changing from -1 at r — to 1 at r — > 00. The first half of (3) is the second order term 
contribution into the mass, the second - the Skyrme term contribution. At fixed r = r the 
latter coincides exactly with 1-dimensional domain wall energy. It is possible to write the 
second order contribution in (2) in the form: 
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and similar for the 4-th order Skyrme term. The equality cj>' = yj2B/A N {\ - 4> 2 )/r eliminates 
considerable part of integrand in (4). Therefore, it is natural to consider function <f> satisfying 
the following differential equation: 

tf = ^-<P) (5) 
which has solution satisfying boundary conditions 0(0) = -1 and 0(oo) = 1: 

(r/r ) b - 1 

with arbitrary r - the distance from the origin of the point where = and profile / = n/2. 
r can be considered as a radius of multiskyrmion. 

After substitution of this ansatz one obtains the soliton mass in the form: 
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Integrating over dr using known expressions for the Euler-type integrals, e.g. 
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allows to obtain the mass of multiskyrmion in simple analytical form as a function of pa- 
rameters b and r : 
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which gives after simple minimization in r 
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and 
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For any value of the power b (10) provides an upper bound for the mass of RM skyrmion. 
To get better bound we should minimize (10) in b. At large enough B when it is possible to 
neglect the influence of slowly varying factors (1 - 1/b 2 ) and bsin(n/b) we obtain easily that 
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and, therefore, 
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The lower bound in (12) was known previously || |(|. The correction to the value bo can be 
found including into minimization procedure the factor (1 - 1/b 2 ) and variation of bsin(n/b) ~ 
tt[1 -tt 2 /(66 2 )]. It provides: 
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and the value b = b + Sb should be inserted into (10). This improves the values of M/B 
for B = 1,2,3... but provides negligible effect for 6—17 and greater, since Sb — 1/Vb. The 
comparison of numerical calculation result and analytical upper bound (12) is presented in 
the Table. 
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The skyrmion mass per unit B-number in universal units 37r 2 _F' 7r /e for RM configurations, approximate and 
precise solutions. The approximate values are calculated using formula (10) with the power b = bo + 8 b. The 
numerical values for 577(2) model are from the papers |3j and earlier papers. The last 3 lines show the result 
for SU(3>) projector ansatz 0] and approximation to this case, An = 4/3 . 



Numerically (10) -(13) provide upper bound which differs from the masses of all known 
RM skyrmions within ~ 2%, beginning with B = 4, see Table. Even for B = 1, where the 
method evidently should not work well, we obtain M = 1.289 instead of precise value M = 1.232. 
For maximal values of B between 17 and 22 where the value of I is calculated, the upper 
bound exceeds the RM value of mass by 0.5% only. We took here the ratio R I/B = 1/B 2 in the 
cases where this ratio is not determined yet, the same as for highest B where it is known, 
i.e. 1.28 for SU(2) caseg, B = 32 and 64, and 1.037 for B > 6 in SU{3) §]. For R I/B = 1 the 
numbers in Table decrease by ~ 0.1% only in the latter case. Note also, that asymptotically 
at large B the ratio of upper and lower bounds 
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~ 1.0396, (14) 



-^7rrm.T 4/2 



M m in 7T\3 y 

i.e. the gap between upper and lower bounds is less than 4%, independently on £?, the 
particular value of 1 and the number of flavors N. With decreasing I the upper bound 
decreases proportionally to the lower bound. In view of such good quantitative agreement 
of analytical and numerical results the studies of basic properties of bubbles of matter made 
in present paper are quite reliable. 

The width (or thickness) W of the bubble shell can be estimated easily. We can define 
the half- width as a distance between points where 4> = ±1/2, then: 
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It is clear that at large B W is constant, and does not depend on the number of flavors N as 
well. The radius of the bubble (11) grows with increasing B like [1/An] 1 ^ 4 . 

Previously we considered large B skyrmion within the " inclined step" approximation 
0. Let W be the width of the step, and r - the radius of the skyrmion where the profile 
/ = ix j2. f = n/2 - (r - r )n/W for r a - W/2 < r < r a + W/2. This approximation describes the 
usual domain wall energy with accuracy ~ 9.5% . 

We wrote the energy in terms of W, r , then minimized it with respect to both of these 
parameters, and find the minimal value of energy. 
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This gave 
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and, after minimization, rl min = y/32/(8A N ) ~ QM2y/l/A N , close to the above result r\ min ~ 
0.667-^/X/Ajv. In dimensional units r = (6X/A A r) 1/4 /(F 7r e). Since X > B 2 , the radius of minimized 
configuration grows as Vb, at least. W min = n, i.e. it does not depend on B for any SU{N), 
similar to previous result (15) which gives W ~ 3.59 for large B, all in units 2/(F„e). The 
energy 

M min ~ (25 + V3AatX/2)/3 (18) 

In difference from previous result (12), (18) does not give the upper bound for the skyrmion 
mass, and for small B, indeed, the value of (18) is smaller than calculated masses of skyrmions. 
For SU(2) model A N = 1 and the energy M min = (2B+y / 31/2)/3. The formula gives the numbers 
for B = 3, 22 in agreement with calculation within RM approximation within 2-3% []3|, |!|. 

More detailed analytical calculation made here confirms the results of such "toy 
model" approximation and both reproduce the picture of RM skyrmions as a two-phase 
object, a spherical bubble with profile / = n inside and / = outside, and a fixed thickness 
shell with fixed surface energy density, p™ rf ~ (2B + y/XTA N /2)/(12m%). The average volume 
density of mass in the shell is 

Pm =-^(2B + ^TA^)^A^JXF^ 2 (19) 

and for SU{2) model at large B it is approximately equal to ~ 0.3Gev/Fm 3 at reasonable choice 
of model parameters F v = 0.186 Gev, e = 4.12 i.e. about twice greater than normal density 
of nuclei. 

3. Consider also the influence of the chiral symmetry breaking mass term which is 
described by the lagrangian density 

M.t. = rh J r 2 (l - cos f)dr, cos f = <f> (20) 

m = 8fi 2 /(3nF 2 e 2 ), n = m w . For strangeness, charm, or bottom the masses m K , m D or m B can 
be inserted instead of m n . 

Instead of above expression (8) we obtain now 

M(B,r ,b) = a(B,b)r + /3(B,b)/r +mrl (21) 

with a, (3 given in (8) and m = 27rm/(6sm(37r/6)). It is possible to obtain precise minimal value 
of the mass 
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When the mass m is small enough, the expansion in 12m/3/a 2 can be made, and one obtains 
the following reduction of dimension r : 
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We used that at large B 



As it was expected from general grounds, dimensions of the soliton decrease with increasing 
m. However, even for large value of m the structure of multiskyrmion at large B remains 
the same: the chiral symmetry broken phase inside of the spherical wall where the main 
contribution to the mass and topological charge is concentrated. The value of the mass 
density inside of the bubble is defined completely by the mass term with 1-0 = 2. The 
baryon number density distribution is quite similar, with only difference that inside the bag 
it equals to zero. It follows from these results that RM approximated multiskyrmions cannot 
model real nuclei at large B, probably B > 12 - 20, and configurations like skyrmion crystals 
H may be more valid for this purpose. 

It is of interest to study what happens at very large values of the mass, when l2m/3 > 
a 2 . Making expansion in a 2 /12mj3 we obtain: 

^-[^(3™)]^, M m ^B,b) ~ ^(3m/3 3 ) 1/4 - (27) 

Minimizing M{B.b) in b gives b^ m = y/zT/B, (3 mm = 8%/2Tb/(9tt) and 
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Af m ^oo(B) ~ y[4V5m/(37r 3 )] 1 / 4 [B2] 3 / 8 . (28) 

So, in this extreme case b — \/B, r — (^)~ 1/2 B 3/8 , the mass of the soliton increases as M — 
^/TXB 9 / 8 , and the binding should become weaker with increasing baryon number of skyrmion. 

It is possible also to calculate analytically the tensors of inertia of multiskyrmion 
configurations within this approximation, see for definitions and formulas. 

4. By means of consideration of a class of functions (6) we established the link between 
toplogical soliton models in rational maps approximation and the soliton models of "domain 
wall" or "spherical bubble" type. Of course, it is a domain wall of special kind, similar to 
honeycomb. The upper bound for the energy of multiskyrmions is obtained which is very 
close to the known energies of RM multiskyrmions, especially at largest B, and is higher 
than the known lower bound by ~ 4% only. 

The following properties of bubbles of matter from RM multiskyrmions are established 
analytically, mostly independent on particular values of the quantity 1: 

The dimensions of the bubble grow with B as \/B, or as X 1/4 , whereas the mass is 
proportional to ~ B at large B. Dimensions of the bubble decrease slightly with increasing 
N - the number of flavors, r - [N/(2(N- 1))] 1/4 , see (11), (12). 

The thickness of the bubbles envelope (15) is constant at large B and does not depend 
on the number of flavors, therefore, the average surface mass density is constant at large 
B, as well as average volume density of the shell. Both densities increase slightly with 
increasing N. At the same time the mass and B-number densities of the whole bubble — > 
when B — ► oo, and this is in contradiction with nuclear physics data. 

It follows from the above consideration that the spherical bubble or bag configuration 
can be obtained from the lagrangian written in terms of chiral degrees of freedom only, i.e. 
the Skyrme model lagrangian leads at large baryon numbers to the formation of spherical 
bubbles of matter, i.e. it provides a field-theoretical realization of the bag-type model. 
This picture of mass distribution in RM multiskyrmions condradicts at first sight to what is 
known about nuclei, however, it emphazises the role of periphery of the nucleus and could be 



an argument in favour of shell- type models of nuclei. The skyrmion crystals || are believed 
to be more adequite for modelling nuclear matter. 

It would be of interest to perform the investigation of the dynamics of bubbles in 
chiral soliton models similar to that performed recently for the simplified two-component 
sigma model, or sine-Gordon model in (3 + 1) dimensions @. Observations concerning the 
structure of large B multiskyrmions made here can be useful in view of possible cosmological 
applications of Skyrme-type models, see e.g. |I0|] . The large scale structure of the mass 
distribution in the Universe Jll[] is similar to that in topological soliton models, and it can 
be the consequence of the similarity of the laws in micro- and macroworld. 

I am indebted to P.M. Sutcliffe for informing me about the results of 0b, c before 
publication and to W.J.Zakrzewski for interest in the questions discussed in present paper. 
The work is supported by RFBR grant 01-02-16615. 
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